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ABSTRACT 

We study the velocities of the accretion along streams from the cosmic web into massive 
galaxies at high redshift with the help of three different suites of AMR hydrodynamical 
cosmological simulations. The results are compared to free-fall velocities and to the 
sound speeds of the hot ambient medium. The sound speed of the hot ambient medium 
is calculated using two different methods to determine the medium’s temperature. 
We find that the simulated cold stream velocities are in violent disagreement with 
the corresponding free-fall profiles. The sound speed is a better albeit not always 
correct description of the cold flows’ velocity. Using these calculations as a first order 
approximation for the gas inflow velocities Unflow = 0.9 i> v ir is given. We conclude from 
the hydrodynamical simulations as our main result that the velocity profiles for the 
cold streams are constant with radius. These constant inflow velocities seem to have a 
“parabola-like” dependency on the host halo mass in units of the virial velocity that 
peaks at M v j r = 10 12 M 0 and we also propose that the best fit functional form for the 
dependency of the inflow velocity on the redshift is a square root power law relation: 

Unflow OC y/Z T I V,j\ r . 

Key words: cosmology: theory - galaxies: evolution - galaxies: formation - galaxies: 
high redshift - intergalactic medium - methods: numerical 


1 INTRODUCTION 


In the recent years it has been shown i n theoretical work 
and i n simulations jFaidal_eta]J 2001 Birnboim_&Dekel 
20031: iKeres et alj 120051. l2009l: iDekel fe BirnboimT 200£; 


Ocvirk. Pichon fe Tevssieil 20081 : Dekel et al. 12009... 

20131') that massive galaxies in the high redshift regime 


(z Xi 1.5), acquire their baryons primarily via cold nar¬ 
row streams of relatively dense and pristine gas with 
temperatures around 10 4 K. These cold flows penetrate 
through the diffuse shock-heated medium. Their activ¬ 
ity peaks around redshift 3. Having reached the inner 
parts of the host halo their gas will eventually form a 
dense, unstable, turbulent disc with a bulge where rap id 
star form ation is t riggered (lAgertz. Tev ssier fe Moorel 


Ceverino. Dekel & Bournaud l2010i: 

iKrumholz & Burkertl 

20ld: Ceverino et al. 

2012 

2014 20151: 

Cacciato. Dekel & Genel 2012 

: Genel. Dekel & Cacciato 


20121 : iGenzel et al.ll201ll : iMandelker et al.ll20l4) . 


* tobias.goerdt@univie.ac.at 


Galaxies’ gas depletion time scales are understood to 
be short wi th respect to the Hubble time throu ghout cos¬ 
mic history dDaddi et al.ll2Qlol : lGenzel et al.ll2010l f. So it be¬ 
comes apparent that galaxies must accrete fresh gas from 
the intergalactic medium to sustain the observed levels of 
star formation over such a long time, since otherwise their 
gas reservoirs would have been completely emptied. In the 
cold flow scenario star formation behaves like a shower that 
is not draining properly: no matter what are the initial con¬ 
ditions, a steady state will be reached at the point where 
the amount of material that comes in via cold streams is 
equal to the drained material, i.e. the gas that is being lost 
by star formation and outflows. In the bathtub it is that 
the pressure at the bottom that rises as it gets more full, 
pushing out more water. For galaxies, it is the star forma¬ 
tion rate that increases with the increase of available gas. 
The result is a self-regulated system, the accretion is driv¬ 
ing the star-formation r ate and not the amount of available 
gas dBouche et al.ll2010l f . It is in this sense that the star for¬ 
mation rate of a galaxy is fundamentally limited by its gas 
inflow and does not depend so much on local physics, such 
as how well the feedback regulates it. The timescales for star 
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2 Tobias Goerdt & Daniel Ceverino 


formation as well as the star formation rate both depend cru¬ 
cially on the gas accretion rate, which in turn depends on 
the speed of the infalling gas when it reaches the halo centre. 
So doubling the infall velocity would also result in doubling 
the star formation rate assuming constant cold stream cross 
sections and densities, which highlights the importance of 
the inflow velocity. 

Processes of observing the cold accr etion stream 
parad igm in the real universe are ongoing. lOoerdt et al.1 
(120121 ') made theoretical predictions about the likelihood of 
observing these streams in absorption. Striking features i n 
absorption have in deed been observed JBouche et al.ll2013l h 
Ilbata et al.l d2013ll found the existence of a planar subgroup 
of satellites in the Androme da galaxy (M 31), compris¬ 
i ng ab out half the population. IGoerdt. Burkert fc Ceverinol 
(120131 ') demonstrated that this vast thin disk of satellites can 
naturally be explained within the cold stream framework 
and should therefore be interpreted as indirect observational 
evidence for the cold stream paradigm. 

lOoerdt et al.l t 201Clh used the cold stream model to ex¬ 
plain the observational phenomenon of Lya-blobs (LABs). 
They deployed cosmological hydrodynamical AMR simula¬ 
tions to predict the characteristics of Lya emission from 
the cold gas streams. These authors pointed out that in 
their simulations the velocities of the cold streams are con¬ 
stant. So the potential energy of the inflowing gas is not 
converted into kinetic energy. Subsequently they concluded 
that it is the release of gravitational energy as the gas is 
flowing down the potential gradient towards the halo cen¬ 
tre which powers the Lya luminosity. The potential en¬ 
ergy of the inflowing gas is released as lya radiation. This 
model was put forward as an alternative to other mod¬ 
els trying to explain the emission of LABs like being pro¬ 
duced in the galaxy’s H 11 regions but then being scattered in 
our direction b y the H 1 gas in the galaxy’s circum-galactic 
medium dSteidel et alJlioilh . be ing lit up due to illumina- 


Scarlata et al.ll200fi 

Prescott et al. 

201511 or being lit up due 

to stellar feedback 

Chapman et al. 

20011: iGeach et al.l[2005l. 


120071') . 


Fauchcr-Gigucrc ct al. (j2010l ~) tested the results of 
IGoerdt et alJ rt201fll ~). They analysed cosmological smoothed 
particle hydrodynamics simulations to predict the Lya cool¬ 
ing emission of forming galaxies wit h the help of a Lya 
radiative transfer code. In contrast to iGoerdt et al.1 ( 2010 1 
they mention that a scenario in which the cold streams free- 
fall into the haloes and release little energy before flitting 
the central disk in a strong shock, may be sup ported by 
their results. Later also lRosdahl fc Blaizotl (120121 ') addressed 
the same question independently. For this purpose they ran 
and analysed cosmological radiation-hydrodynamics simu¬ 
lations. They looked at phase diagrams for their haloes of 
gas speed versus radius in which they were able to distin¬ 
guish several cold streams. They produced corresponding 
free-fall profiles (i.e. the speeds that the strea ms would fol¬ 
low if they were in free-fall). In contrast to iGoerdt et al.l 
(I2010T) they concluded that the streams are close to free- 
fall. As we will also discuss in greater depth in section 0 
the discrepancies between their work and the results pre- 
sented later in this paper m ight result from the fact, that 
iFaucher-Giguere et al.l J2010I *) on the one hand deployed SPll 
simulations which used to have inaccuracies modelling dis¬ 


sipative p rocesses (|Agertz__eT_alJ |2007|) and the sample of 
galaxies of lRosdahl fe Blaizotl ( 20121 1 on the other hand was 
very small (only three galaxies in total) and they report that 
not all of their galaxies showed signs of free-fall behaviour. 

Rad ial velocities around 10 1 2 haloes have been 
studied (Ivan de Voort fc Schavdl2012fl . They found roughly 
constant values of the radial inflow velocities, with a decrease 
tow ards smaller r adii. 

iNelson et al.l (120131 1 compare results from SPH with 
moving mesh simulations and find that the filamentary ge¬ 
ometry of accreting gas near the virial radius is a common 
feature in massive haloes above 10 11,5 Mg. Gas filaments in 
GADGET, however, tend to remain collimated and flow coher¬ 
ently to small radii, or artificially fragment and form a large 
number of purely numerical blobs. These same filamentary 
gas streams in AREPO show increased heating and disruption 
at 0.25 - 0.5 r v i r and contribute to the hot gas accretion rate 
in a manner distinct from classical cooling flows. 

INelson et al.1 (l2015al l investigated how the way galaxies 
acquire their gas across cosmic time in cosmological hydro- 
dynamic simulations is modified by a comprehensive physi¬ 
cal model for baryonic feedback processes finding that feed¬ 
back strongly suppresses the raw, as well as the net, inflow 
of this smooth mode gas at all redshifts, regardless of the 
tem perature history o f newl y acquired gas. 

I Wetzel fc Nagail (120141 ') explored systematically the 
physics that governs cosmic accretion into haloes and their 
galaxies. They analysed a suite of cosmological simulations 
that incorporate both dark matter and gas dynamics with 
differing treatments of gas cooling, star formation and ther¬ 
mal feedback to explore the physics that governs the accre¬ 
tion of dark matter and baryons into halos and their galax¬ 
ies. They presented average radial velocities of the inflowing 
material (their figures 4 and 6, bottom panels). The inflow 
velocities presented there decrease with decreasing radius 
(roughly a decrease of 65% between 1 r v n and the centre of 
the host)._ 

In a companion paper (IGoerdt et al.ll2015l f we look at 
the amount of inflow - the mass accretion rate - both as 
a function of radius, mass and redshift for the three con¬ 
stituents gas, stars and dark matter. 

In this paper we look at the velocity profiles of the 
gas inflow, both as a function of radius, mass and redshift 
for the three constituents of the simulations, namely gas, 
stars and dark matter. This study is worthwhile doing since 
multiphase gas accretion is an interesting topic in galaxy 
formation theory by itself and this study in particular has 
important implications on the judgement of plausibility of 
both star formation rates as well as the gravitational cooling 
model for LABs, as already mentioned. The paper is organ¬ 
ised as follows: in section [2] we present the three suites of 
simulations used for the following analysis. In section [3] we 
show a brief analysis of the mass accretion rate of the gas. In 
section 0 we present our toy models. In sectior0we discuss 
the velocity profiles of the inflow along the streams and in 
section 0 we draw our conclusions. 


2 SIMULATIONS 

We use snapshots of galaxies from three different sets of sim¬ 
ulations, all three employing Eulerian AMR (Adaptive Mesh 
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Inflow velocities of cold streams 3 



Figure 1. Two simulated galaxies from ART and MN. The colour refers to the maximum gas density along the line of sight. Contour 
lines are at n = 0.1, 0.01 and 0.001cm -3 , respectively. The circles show the virial radii. Left: a typical ART galaxy (resolution 70pc) at 
z = 2.3, with M v i r = 3.5 X 10 11 Mg. Right: one of the MN galaxies (resolution 1 kpc) at z = 2.5, with M v j r = 10 12 Mq. In both cases, 
the inflow is dominated by three cold narrow streams that are partly clumpy. The density in the streams is n = 0.003 — 0.1 cm -3 , with 
the clump cores reaching n ~ 1cm -3 . 


Refine ment) hydrodynamics in a cosmolog ic al setting. The 


ART |Ceverino. Dckel & Bournaud 20101: 

Ceverino et al.l 

20121: Dckel et al. 20131; Oeverino et al.1 2015 

) suite consists 


of several zoom-in simulations with a maximum resolution 
of 35 — 70 pc at z = 2. The simulation zooms in on in¬ 
dividual galaxies that reside in dark-matter haloes which 
have masses of (0.13 — 1.60) x 10 12 Mq at z = 2.3. This 
ART suite has got an extension that we will also use, the 
Arp suite of simulations which includes radiation pressure 
feedback (ICeverino et al.ll2014l ). This second suites consist 
of individual galaxies that reside in dark-matter haloes of 
masses (1 — 8) x 10 11 Mq at z = 2.0. The third simula¬ 
tion is the Horizon-MareNostrum (h ereafter MN) simula¬ 
tion (lOc.virk. Pichon fc Tevssierll2008i ) containing dozens of 
massive galaxies in a cosmological box of side 71 Mpc with 
a maximum resolution of ^ 1 kpc . 

Density maps of galaxies of the MN as well as the ART 
suite are shown in figurc[T| They demonstrate the dominance 
of typically three, narrow cold streams, which come from 
well outside the virial radius along the dark-matter filaments 
of the cosmic web, and penetrate into the discs at the halo 
centres. The streams are partly clumpy and partly smooth, 
even in the simulation of higher resolution. The typical den¬ 
sities in the streams are in the range n = 0.01 — 0.1cm -3 , 
and they reach n = 0.1 — 1cm -3 at the clump centres and 
in the central disk. 


2.1 High-resolution ART simulations 

These simulations wer e run with the AMR code ART (Adap¬ 
tive Refinemen t Tree: iKravtsov . Kly pin fc Kho khlov! 1 19971 : 
lKravtsovll2003l ) with a spatial resolution better than 70 pc in 
physical units. It incorporates gas cooling, photoionisation 


heati ng, star formation, metal enrichment and stellar feed¬ 
back dCeverino fe Klypinll2009l ) which are relevant physical 
processes for galaxy formation. For the gas density, tem¬ 
perature, metallicity, and UV background in any given cell 
coolin g rates were computed using CLOUDY (IFerland et al.l 
1 1998h . It is assumed that coolin g takes place at the cen¬ 
tre of a cloud of thickness 1 kpc (ICeverino-Rodrigueal2008l : 
iRobertson fe Kravtsovll2008l ) . Metallicity dependent, metal¬ 
line cooling is also included, assuming a relative abundance 
of elements equal to the solar composition. The code im¬ 
plements a “constant” feedback model, in which the com¬ 
bined energy from stellar winds and supernova explosions 
is released as a constant heating rate over 40 Myr (the 
typical age of the lightest star that can still explode in 
a type-II supernova). Photo-heating is taken into account 
self-cons istently with radia t ive co oling. A uniform UV back¬ 
ground (lHaardt fe Madaul 1 19961 ) is assumed, ignoring lo¬ 
cal sources. The self-shielding of dense, galactic neutral 
hydrogen from the cosmological UV background is mim¬ 
icked by assuming a substantially suppressed UV back¬ 
ground (5.9 x 10 26 erg s _1 cm -2 Hz -1 , the value of the 
pre-reionisation UV background at z = 8, according to 
the Haard & Madau model) for the gas at densities above 
n = 0.1 cm -3 . This density threshold is mo tivated by radia¬ 
tive transfer calculations (appendix A1 of iFumagalli et all 
l201lL figure A2). 

Our special version of the ART code has a unique fea¬ 
ture for the purpose of simulating the detailed structure of 
the streams. It allows gas cooling to well below 10 4 K so 
enabling the formation of high densities in pressure equi¬ 
librium with the hotter and more dilute medium. A non- 
thermal pressure floor has been implemented to ensure that 
the Jeans length is resolved by at least seven resolution ele- 
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men 

est g 

s and thus prevent artificial fragmentation on the small- 

rid scale dTruelove et al. 1997 

: Robertson & Kravtsovl 

2008 

Ceverino. Dckcl & Bournauc 

2010h. It is effective in 


the dense (n > 10 cm 3 ) and cold (T < 10 4 R) regions in¬ 
side galactic disks. 


The equation of state remains unchanged at all den¬ 
sities. Stars form accordin g to a stochastic model t hat is 
described in appendix A of I Ceverino fe Kiypinl d2009l ). The 
SFR. surface density and the gas surface density are mea¬ 
sured from the simulations and co mpared with the Kenn i- 
cutt relation in the appendix A of I Ceverino et ahl J2014ll . 
figu re Al . The s tocha stic star formation model reproduces 
the iKennicu tt] Jl998l l law. This happens in cells where 
the gas temperature is below 10 4 K and the gas density is 
above a threshold of n = 1cm -3 . We use the IMF from 
Miller fc Scald (Il979ll tha t is also matching the results of 
Wooslev fc Weaver! (| 19951 1. Metals from supernovae type II 
and type la enrich the ISM. They are released from each 
star particle by SNII at a constant rate for 40 Myr after its 
birth. The metal ejection by SNIa assumes an exponentially 
declining SNIa rate from a maximum at 1 Gyr. The code 
treats the advection of metals self-consi stently and it distin¬ 
guish es between SNII and SNIa ejecta (ICeverino-Rodriguea 

l20oah . 

2.1.1 Simulations including radiation pressure (Aycp) 

The Arp suite of simulations (ICeverino et al.l l2014l l is a 
further development of last subsection’s ART suite: apart 
from the features already presented there, it also includes 
the effects of radiation pressure by massive stars. The ra¬ 
diation pressure was modelled as a non-thermal pressure 
that acts only in dense and optically thick star-forming re¬ 
gions in a way that the ionising radiation injects momen- 
tum around massiv e stars, pressurising star-forming regions 
jAaertz et al.ll2013l . their appendix B). The adaptive comov¬ 
ing mesh has been refined in the dense regions to cells of 
minimum size between 17-35 pc in physical units. The DM 
particle mass is 8.3 x 10 4 Mq. The particles representing 
star clusters have a minimum mass of 10 3 Mq, similar to 
the stellar mass of an Orion-like star cluster. 


2.2 Ramses Horizon-MareNostrum simulation 

The AMR code RAMSES dTevssierl l2002ll is used for this 
simulation. Its spatial resolution in physical units is 
~ 1 k pc. We included UV heating via a background 

model dHaardt &: Madaul 1 1996l l. Supernovae feedback and 
metal enrichment are modelle d in a simple wav usin g 
the implementation described in lDubois fc Tevssier <2008 f). 
Cooling rates are calculated assuming ionisation equilib- 
riurn for H and He, for coll i siona l- and photo-ionisation 
dKatz, Hernauist fe Weinberd Il992l l. Metal cooling is as¬ 
sumed to be pro portional to the metallicit y, relative to the 
solar abundance dGrevesse fo Sauvalll998l l. It is calculated 
with the help of tabulated CLOUDY rates. Unlike in the ART 
or Arp simulations, no cooling below T < 10 4 K is computed, 
and no self-shielding of the UV flux is assumed. 

For high-density regions, we consider a polytropic 
equation of state with 70 = 5/3 to model the com¬ 
plex, turbulent multi-phase structure of the inter-stellar 


label 

suite 

M vir [10 12 Mq] 

Z 

Ngal 

10 13 

MN 

10.47 ± 0.56 

1.57 

12 

10 13 

MN 

10.49 ± 0.93 

2.46 

12 

5 x 10 12 

MN 

5.00 ± 0.045 

1.57 

12 

5 x 10 12 

MN 

5.48 ± 0.26 

2.46 

11 

10 12 

MN 

1.03 ± 0.003 

1.57 

8 

10 12 

MN 

1.01 ± 0.004 

2.46 

12 

10 12 

MN 

1.03 ± 0.006 

4.01 

9 

10 11 

MN 

0.099 ± 0.000 

1.57 

12 

10 11 

MN 

0.099 ± 0.000 

2.46 

7 

10 11 

MN 

0.099 ± 0.000 

4.01 

12 

1.9 x 10 12 

ART 

1.907 ± 0.217 

1.14 ± 0.02 

34 

1.3 x 10 12 

ART 

1.286 ± 0.093 

1.60 ± 0.02 

73 

8.6 x 10 11 

ART 

0.863 ± 0.046 

2.25 ± 0.02 

109 

3.9 x 10 11 

ART 

0.391 ± 0.034 

3.40 ± 0.04 

119 

7.1 x 10 11 

Arp 

0.707 ± 0.055 

1.14 ± 0.02 

41 

6.7 x 10 11 

Arp 

0.672 ± 0.049 

1.58 ± 0.02 

47 

4.9 x 10 11 

Arp 

0.491 ± 0.028 

2.21 ± 0.03 

57 

2.9 x 10 11 

Arp 

0.290 ± 0.019 

3.27 ± 0.05 

62 

2.6 x 10 11 

Arp 

0.260 ± 0.010 

1.13 ± 0.02 

40 

2.4 x 10 11 

Arp 

0.242 ± 0.010 

1.59 ± 0.02 

50 

1.6 x 10 11 

Arp 

0.160 ± 0.007 

2.29 ± 0.03 

60 

7.3 x 10 10 

Arp 

0.073 ± 0.004 

3.53 ± 0.05 

66 


Table 1. The different bins of galaxies which are used through¬ 
out this paper. ’Label’ denotes the tag a bin is labelled with in 
any of the figures. It usually is a mass very close to the ensem¬ 
ble’s actual mean virial mass. ’Suite’ indicates the simulation set 
the bin consists of. M v ; r gives the actual mean virial mass of the 
ensemble together with its standard deviation, z is the mean red- 
shift of the ensemble together with its standard deviation. We do 
not combine galaxies from different MN snapshots therefore the 
standard deviation of the redshift of any of the MN bins is always 
zero. lVg a i is the number of galaxies in the respective bin. 


medi um (ISM) CYepes et al. 1997; ISnringel fc Hc rnciuistl 
2003|) in a simplifie d form (see Schave fc Dalla Vecchiall2008l : 
Dubois fe Tevssierll2008l 'l. The ISM is defined as gas with hy¬ 
drogen density greater than uh = 0.1cm -3 . This definition 
is one order of magnitude lower than the one of the ART 
and Arp simulations. Star formation only for ISM gas has 
been included. It is mod elled by spawning star particles at a 
rate consistent with the iKennicuttl (Il998ll law derived from 
local observations of star forming galaxies. 

The RAMSES code implements a pressure floor in or¬ 
der to prevent artificial fragmentation, by keeping the Jeans 
lengthscale, Aj oc Tn ~ 2 ^ 3 , larger than four grid-cell sides 
everywhere. In any case where n > 0.1cm -3 , a density de¬ 
pendent temperature floor was imposed. It mimics the aver- 
age thermal and turbulent pressure of the multip ha se ISM 
dSpringel fe Hernauistll2003l ; iDalla Vecchia fe Schavdl^OOSl l . 

In our case, we allow the gas to heat up above this tem¬ 
perature floor and cool back. The temperature floor follows 
a polytropic equation of state with Tfl oor = To(n/no) 70-1 , 
where To = 10 4 K and no = 0.1 atoms cm -3 . The resulting 
pressure floor is given by Pfl oor = nHfeIfl oor . For each stel¬ 
lar population, 10% of the mass is assumed to go through 
a supernovae type II event after 10 Myr. For each super¬ 
nova, 10% of the ejected mass is assumed to be pure metals. 
The energy and the metals produced by the supernova event 
are released in a single impulse. The remaining 90% of the 
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Figure 2. Average inflow as a function of radius for different halo masses and redshifts for the Arp simulations. The different panels 
indicate a variety of different halo masses and redshifts. The inflow in gas is remarkably constant with radius, except for the higher 
mass very low redshift panel, which has a very mildly declining mass accretion with increasing radius. However by and large the mass 
accretion rate of gas is constant with radius. 


ejected mass keep the metallicity of their star at birth. Feed¬ 
back through supernovae type la has not been considered. 

The initial conditions were constructed assuming a 
ACDM universe with 12m = 0.3, S2a = 0.7, fib = 0.045, 
h = 0.7 and as = 0.9 in a periodic box of 71 Mpc. 
The adaptive-resolution rules in this simulation were the 
same everywhere, with no zoom-in resimulation of individ¬ 
ual galaxies. The dark matter particle mass is 1.16 x 10 7 Mq, 
the star particle mass is 2.05 x 10 6 Mq, the smallest cell size 
is 1.09 kpc physical, and the force softening length is 1.65 
kpc. 

We will show averaged results for an ensemble of galax¬ 
ies having very similar masses and redshifts. In table Q] we 
show a summary of the various bins of galaxies we use. Since 
we have better statistics for the MN simulation we could bin 
galaxies from a narrower mass range therefore the standard 
deviations of the mean mass of a MN bin is usually much 
smaller. To partly compensate for that we combine ART 
as well as Arp galaxies from adjacent redshifts increasing 
the statistics but introducing a standard deviation into the 
mean redshift of the bin. 


3 MASS ACCRETION RATE 

In a companion paper dGoerdt et al.|[2015l l we look at the 
amount of inflow - the mass accretion rate - both as a func¬ 
tion of radius, mass and redshift for the three constituents 
gas, stars and dark matter. Since the mass accretion rate of 
the gas is ultimately what is interesting for the rate of star 
formation and how much gas needs to be lost via outflows 


we show briefly an analysis of the mass accretion rate of the 
gas in this section. 

The analysis of the Arp simulations starts by looking at 
the average gas inflow as a function of radius: We measure 
the amount of mass which is crossing a spherical shell of 
radius r centred around a given individual galaxy within a 
small time At. The mass crossing the shell is divided by the 
time taken (At) to get an inflow rate. In order to get the best 
estimate for the inflow in case of gas we decided to take into 
account only the cells which have an inwards radial velocity. 
In order to get rid of the statistical noise of a single galaxies 
we average the amount of inflow over all available galaxies 
having similar redshifts and masses. 

The amount of inflow into the galaxies as a function 
of radius is shown in figure [2] The magnitude of the inflow 
is basically constant over all radii r regardless of the ac¬ 
tual mass-redshift bin, except for the one higher mass very 
low redshift panel, which has a very mildly declining mass 
accretion with increasing radius. However by and large the 
amount of inflow is constant with radius. So we can mainly 
see that the mass accretion rate simply follows the velocity 
profile, i.e. cold and hot gas are not mixed as it falls in. As 
already mentioned we discuss these trends and also the mass 
accr etion rates of stars and dark matter in a companion pa¬ 
per dGoerdt et al.ll2015h . 


4 ANALYTICAL MODELS 

Naively one would expect the gas inflow velocities to obey 
one of the following two quantities: either free-fall or the 
sound speed of the hot ambient medium. The underlying as- 
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sumption of the first case is that the gas flows freely within 
the gravitational potential into its host and is otherwise com¬ 
pletely unconstrained. In the second case one assumes that 
the velocity of the gas will be fundamentally limited by the 
sound speed of the hot ambient medium through the effect 
of shocks which occur within the streams themselves as well 
as at the boundary layers between the streams and the hot 
ambient medium as soon as the velocity of the inflowing gas 
reaches the sound speed of this ambient medium. 


4.1 Free-fall 


Free-fall profiles cannot be calculated analytically. To eval- 
uate them numeri c ally w e follow the method laid out by 
iRosdahl fc Blaizotl (I2012T ) which we will quickly outline here: 
a static, spherically symmetric potential is assumed. We in¬ 
tegrate the free-fall speeds vs(r) from a starting position ro 
and a starting velocity Uinit = Vfr(ro) towards the halo centre 
using 


dv ff(r) 


1 GM « r ) dr , 

v a (r) r 2 


(1) 


where r is radius, M{< r ) is the total halo mass within r 
and G is the gravitational constant. In practice we divide the 
host halo into spherical bins at radial positions n, where an 
increasing i corresponds to a decreasing radius, and solve 
the above equation by recursively computing 


Vff(n+i) 


vff(n) + 


1 G M{< n+i) 

/ \ 2 v i+1 

v{n) rf +1 


n). 


( 2 ) 


The underlying mass distribution M(< n+i) is measured 
from the outputs of the hydrodynamical simulations as de¬ 
scribed in section [2] We choose an initial position (ro = 
1.7 r v ir) as well as an initial inward velocity (limit = 0) and 
integrate the path of the free-falling gas iteratively until a 
very small radius (0.1 r v i r ) is reached. The free-fall profiles 
usually show a strongly increasing velocity with decreasing 
radius. The results of the calculations are presented in detail 
in section O We will see that these predictions violently dis¬ 
agree with the results from the hydrodynamical simulations. 


4.2 Sound speed of the ambient medium 


The sound speed c B of the hot ambient halo gas on the other 
hand can be calculated analytically. Sound speed as such 
equals to 


Cs (r) 


7 A:b T(f) 

TTlp av 


(3) 


with 7 being the ratio of the heat capacity at constant pres¬ 
sure to the heat capacity at constant volume, which is usu¬ 
ally a constant (~ 5/3). ks is Boltzmann’s constant, m p 
is the proton mass, £i av is the mean molecular weight (for 
ionised gas /r av = 0.61) and T(r) is temperature of the gas 
as a function of position. 


4-2.1 Approximated constant halo gas temperature 

For a first order approximation the temperature of the ambi¬ 
ent halo gas is assumed to be equal to the virial temperature 
of the host halo. It is constant with position and given by 


rp _ G Alyir m p //;>.v , . 

vir “ 2 r vir k B ' ( } 

G is the gravitational constant, M v i r is the virial mass of the 
host halo and r v ; r its virial radius. These three quantities are 
related via 


Uvir 




G M v ii 

rvir 


(5) 


After inserting the expressions for T v i r and « v ir into equation 
© we get an expression for the sound speed. Assuming that 
the velocity of the gas in the cold streams indeed equals to 
the sound speed of the ambient medium, because the shocks 
which occur as soon as the inflow velocity reaches the sound 
speed will fundamentally limit the velocity of the gas one 
gets a first order approximation for the infall velocity: 


i-’inflow — 0.913 Uvir 


( 6 ) 


In this approximation the infall velocity is constant with 
radial distance from the host halo. We will come back to 
this statement at the end of section [5] where we will see 
that this particular equation tends to be accurate for very 
high redshift haloes (z X! 4.0) and more importantly that the 
infall velocity is in general indeed constant with radius. 


4-2.2 Simulated halo gas temperature 

A more sophisticated measure of the ambient medium’s tem¬ 
perature is not to approximate the gas temperature any¬ 
more, but to measure it directly from the hydrodynamic 
simulations as a function of position. Since we assume spher¬ 
ical symmetry it is sufficient to measure it as a function of 
radial distance from the host halo centre. When looking at 
the simulations one has to carefully discriminate between 
the gas of the ambient medium which is hot and the gas 
flowing into the host as part of the cold streams which is 
cold. Only the temperatures of the ambient medium will 
determine the velocity of the cold streams. To distinguish 
between cold streams on the one hand and hot halo gas 
on the other hand we adopt the following temperature cut: 
temperatures are defined as ’hot’ in this sense if they are at 
least > 0.2 T v i r or > 3 x 10 4 K whatever is the higher value. 
The mass weighted average temperatures of the hot ambient 
halo gas only are then taken in spherically symmetric bins 
at different radii r to get a numerical temperature profile 
T(r). Inserting this numerical profile for the temperature 
into equation © we get a sound speed profile for each set 
of simulated galaxies that are listed in table |T| The sound 
speed profiles are presented together with the results from 
the hydrodynamical simulations in section 0 For the MN 
simulations at intermediate to high masses (M v i r > 10 12 
Mg), the sound speed is increasing with decreasing redshift 
almost as sharply as the free-fall profiles. In these regions 
of the parameter space it can only serve as a slightly better 
proxy for the inflow velocity than free-fall. In all four bins 
of the ART simulations the sound speed profiles increase 
slightly with decreasing radius. In case of the ART simula¬ 
tions the sound speeds match the simulations very precisely 
as we will see later in figure [4] The profiles increase slightly 
with decreasing radius also in the Arp simulations. Their 
shape is very similar to the actual simulation values but the 
sound speeds are usually 50% higher than the simulation 
values. To summarise: the sound speed of the hot ambient 
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Figure 3. Mass weighted averages of the inflow velocities as a function of radius for the MN simulations. The different panels indicate 
a variety of different halo masses and redshifts. Corresponding free-fall profiles are over-plotted in dotted magenta and profiles of the 
sound speed of the hot ambient halo gas in dash-dotted turquoise. For low mass haloes (M v 1T ^ 10 12 Mq) at all redshifts the velocities 
are roughly constant with radius. For high mass haloes (M v i r ) 5 x 10 12 Mq) the inflow velocities increase slightly with increasing 
radius. Subsequently no sign of free-fall can be seen in any of the bins. In the very high redshift and low mass cases one sees an increase 
of velocity with decreasing radius down to 0.6 r v i r from where the velocity is decreasing again. These shapes are typical for high redshift 
low mass haloes. In units of the virial velocity the average inflow velocity stays roughly constant with varying host halo mass, it increases 
with an increase in redshift. 
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Figure 4. Mass weighted averages of the inflow velocities as a function of radius for the ART simulations. Free-fall profiles are over¬ 
plotted in dotted magenta and the sound speeds of the hot ambient halo gas in dash-dotted turquoise. The inflow velocities of the gas 
increase slightly with decreasing radius. The total increase over the whole radius range is usually roughly 50%. These increases however 
are not consistent with free-fall as the corresponding profiles show. The sound speeds of the hot halo gas are a much better proxy, at 
least at high redshift. In the higher redshift cases one sees an increase of velocity with decreasing radius down to ~ 0.5 r v ; r from where 
the velocity is decreasing again. A shape like that is typical for lo w mass haloes at high redshift, it originates in the loss of angular 
momentum during the “strong-torque phase” IIDanovich et al.ll2015l l. 


halo gas is a better proxy for the velocity of the infall than 
free-fall. It is more precise at low masses (M v n ^ 10 12 Mq) 
and gives very good results for the ART suite of simulations. 
However it is only an incomplete prescription for the inflow 
velocity. 


5 INFLOW VELOCITIES 

The analysis of the hydrodynamical simulation starts by 
compiling averaged infall velocity profiles for all simulated 
galaxies listed in table [T] Since we want to look at the ve¬ 
locities of the inflowing cold streams and not at the average 
velocity of some net inflow we have to discriminate between 
the inflowing material belonging to the cold streams and 
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Figure 5. Mass weighted averages of the inflow velocities as a function of radius for the Arp simulations. Free-fall profiles are over¬ 
plotted in dotted magenta and the sound speeds of the hot ambient halo gas in dash-dotted turquoise. Again the inflow velocities of 
the gas increase with decreasing radius in most cases. These increases are again not consistent with free-fall. In most cases one sees an 
increase of velocity with decreasing radius down to ~ 0.5 r v ; r from where the velocity is decreasing again. A shape like that is typical 
for low mass haloes at high redshift. 


the ambient halo material. We experimented with different 
cuts and found out that a simple inflow criterion should be 
sufficient. Since different gas phas es are reported to have 
very different accretion properties dNelson et al.lfioi^ 'l and 
we want to concentrate on the kinematics of the cold flows 
only we decided to adopt for the gas an additional temper¬ 
ature cut of T < 3 x 10 4 K. So the velocities are computed 
for gas by averaging the radial velocity component of all in¬ 
flowing cells having temperatures less than 3 x 10 4 K in a 
given radius bin. The dark matter and stellar velocities are 
computed for by averaging the radial velocity component of 
all inflowing particles in a given radius bin. Only cells and 
particles with negative radial velocities and only cells with 
temperatures less than 3 x 10 4 K are taken into account, 
i.e. all motion perpendicular to the radial infall direction as 
well as all cells or particles which are moving outwards as 
well as all cells with temperatures higher than 3 x 10 4 are 
completely neglected. The resulting values are averaged over 
all available galaxies in the respective mass and redshift bin. 
Velocities are always quoted dimensionless, namely in units 
of the virial velocity u v ; r which has already been defined in 
equation ©• 

The resulting inflow velocity profiles are shown in fig¬ 
ures Q5] (MN), |4] (ART) and [5] (Arp). Free-fall profiles are 
over-plotted in dotted magenta and the sound speeds of the 
ambient gas is over-plotted in dash-dotted turquoise. The 
most striking feature visible is that the velocity is to a very 
good approximation constant with radius. It is clearly not 
free-falling. This is the case for all three suites of simula¬ 
tions, for all redshifts or masses as well as for all three con¬ 
stituents: gas, stars or dark matter. The potential energy of 


the inflowing gas is not converted into kinetic energy, it has 
to be dissi pated by another m echanism, most likely by lya- 
radiation (iGoerdt et alj|201(ih . In some of the cases however 
the velocity of the gas follows the sound speed profiles of the 
hot ambient medium fairly well. 

Looking at the graphs more carefully the following devi¬ 
ations from constant velocity infall can be found: MN’s high 
mass haloes (M v r K x 10 12 Mq) have decreasing inflow 
velocities with decreasing radius. In most of the ART and 
Arp galaxies on the other hand an increase of the velocity 
can be seen as the material falls further in. This infall how¬ 
ever is far from being consistent with free-fall as the overlain 
profiles demonstrate. In all other cases the inflow velocity is 
indeed constant. 

Additionally there is a dependency of the shape of the 
infall velocity profiles with mass and redshift: For low mass 
haloes (M v ir ^ 5 x 10 12 Mg) at high redshift (z Xi 1.7) usu¬ 
ally there is an increase of velocity with decreasing radius 
down to ~ 0.5 r v i r from where the velocity is decreasing 
again. This behaviour originates in the loss of angular mo¬ 
mentum due to strong torques in a til ted extended ring, the 
“strong-torque phase” of the inflow (IDanovich et al.l 20151 . 
their section 5 and figure 20). For very high mass haloes 
(Mvir > 2.5 x 10 12 Mq) the inflow velocity is slightly in¬ 
creasing with increasing radius over the whole radius range 
without ever turning down again. In the remaining bins the 
inflow velocity is constant or only very mildly increasing 
with decreasing radius also over the whole radius range also 
without ever turning down again. 

The velocities of the star particles with negative radial 
velocities stay remarkably constant over the whole radius 
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Figure 6. Direct comparison between the MN (dotted blue), the 
ART (solid red) and the Arp (long-dashed green) suite of sim¬ 
ulations. The inflow velocities as a function of radius are shown 
together with their 1 a standard deviation lines. Over-plotted are 
a free-fall profile (short-dashed magenta) and the sound speed of 
the hot ambient halo gas c s (dash-dotted turquoise). Inflow veloc¬ 
ities for MN, ART and Arp roughly agree with each other within 
the 1 a standard deviation, small offsets can be accounted for by 
the slightly different redshifts. 


range for all masses and redshifts. The velocities of the dark 
matter particles with negative radial velocities increase with 
increasing radius for high mass galaxies (> 3 x 10 12 Mq) 
and decrease with increasing radius for low mass galaxies 
(< 3 x 10 12 Mq). One would naively expect that at least 
the dark matter particles inflow velocities would follow free- 
fall, since those are collision-less and therefore cannot be 
slowed down by any baryonic processes. However there is 
no meaningful dark matter accret i on w ithi n ~ 2 r v j r as al¬ 
ready pointed out bv ICuesta et al.l ll2008l) or lWetzel fe Nagail 
(120141) . Therefore the dark matter particles are not inflow¬ 
ing within the radius range analysed, they are just moving in 
and out again of a certain innermost sphere since the halo is 
virialised and most of the dark matter particles are orbiting 
on non-circular orbits. 

The finding of const ant velocity inflow is_in 

contrast to earlier wo rks dFaucher-Giguere et alj l20ld : 
iRosdahl fe Blaizotl 120121 ) . These authors report evidence 
for free-fall behaviour of the cold s treams. A possible ex¬ 
planat ion for the discrepancies with iFaucher-Giguere et al.l 
ll2010l) might be, that they deployed SPH simulations which 
used to have inacc uracies modelling dissipative processes 
(lAgertz et al.l 120071 ). Such processes are taking place fre- 
quently at t he b oundary layers of cold streams. Also 
I Nelson et al.l d2013l ) found differe nces in the accretion prop - 
erties for different hydro schemes. [Rosdahl fe Blaizotl (120121) 
on the other hand used AMR techniques, like us. Their sample 
of galaxies however was small (only three galaxies in total) 
and they report that not all of their galaxies showed signs 
of free-fall behaviour. Constant velocity inflo w on the othe r 
hand is in agreement with the results from iTeklul d2012l) . 
She finds constant velocity profiles at around 170kms _1 for 
three galaxies with M v i r = 1.5 x 10 12 Mq at z = 2.33 in 
her simulations. Roughly co nstant velocity inflow is also i n 
agreement with the results of ivan de Voort fc Schave] 1201 2). 
their absolute velocities of 150 km/s (~ 0.75 v v ir) is re- 
markably clo s e to th e values we found (~ 0.7 u v ir). Also 
iNelson et al.l ( 2015b ) find roughly constant inflow veloci¬ 
ties, c.f their figure 7. Strictly non f ree-falling inflow i s also 
in agreement with the results from I Wetzel fe Nagail (120141 . 
their figures 4 and 6) who see a decrease in velocity with 
decreasing radius (roughly a decrease of 65% between 1 r v i r 
and the centre of the host). Unlike us they did not exclude 
outflowing cells and particles which explains their velocities 
that are varying with radius as well as the on average lower 
values of their speeds. 

To guarantee the consistency of the MN, the ART and 
the Arp simulations it is necessary to compare them di¬ 
rectly. In figure [6] galaxies at z ~ 2.3 with M v i r ~ 8 x 10 11 
Mq are shown for the three suites of simulations. The lines 
for stars and dark matter are very similar for all three: the 
velocities are roughly constant over the whole radius range 
at comparable values. However the gas curves are slightly 
different: the MN one is by and large constant, whereas the 
ART as well as the Arp one increase slightly with decreasing 
radius. The overall increase is roughly 35 % over the whole 
radius range (1.7 — 0.2 r v r). These increases, as mentioned 
earlier, are not consistent with free-fall. There is an offset in 
the values: the MN velocities are ~ 30% higher over most 
of the radius range. The offset can be accounted for by the 
offs ets, which the host h aloes’ masses and redshifts have. 

INelson et al.l (l2015al ) found that feedback affects gas ac- 
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Figure 7. Average inflow velocity r>i n fl ow in units of the virial velocity as a function of virial mass M v ; r (left) or as a function of redshift 
2 + 1 (right). There is a strong square root power law relationship between the inflow velocity and the redshift (right panel). The inflow 
velocity also follows a “parabola-like” relation with respect to the host halo mass on the other hand that peaks roughly at M v ; r = 10 12 
Mg (left panel). The data points within each panel were scaled with the help of the scaling relation in the respective other panel to the 
values indicated. Both scaling relations were found together in an iterative process (see text). The colour bar axes indicate the values 
the data points used to have before rescaling. 


cretion. Our three suites of simulations use varying strengths 
of feedback. As we can see from figure [6] varying strengths of 
feedback affect the inflow velocities only very mildly. How- 
ever as we w ill discu ss in our forthcoming companion paper 
(iGoerdt et al.|[2015l ) varying strengths of feedback do affect 
the total amount of inflow. 

After the consistency of the infall velocities in all three 
suites of simulations has been established one can now look 
at possible trends of the infall velocities with mass and red¬ 
shift. We want to quantify the qualitative statements made 
earlier. For simplicity a constant infall velocity with radius 
is assumed. Its actual value is determined by a constant line 
fit through the velocity curve of the respective mass-redshift 
bin. 

The average inflow velocities Vmflow of the gas are shown 
as a function of halo mass M v ; r and redshift 2 +1 in figure [7] 
Plotted are in both panels all available bins of galaxies from 
MN, ART and Arp as listed in table [T] In the left panel the 
^inflow values are plotted against the halo mass (M v i r ) and 
in the right panel against redshift 2 +1. One can easily iden¬ 
tify that the Vmfiow values in the right panel follow a linear 
relation with respect to z + 1 on this log-log graph, which is 
a power law. Since in this case the exponent is remarkably 
close to a half we decided in order to reduce the number of 
free parameters of our fitting function to choose a square 
root power law as the fitting function for the Vi n fl ow values 
with respect to 2 + 1. The best fit square root power law is 


over-plotted by the solid black line. The ^inflow values in the 
left panel however seem to follow a “parabola-like” function 
on this log-log graph, which is the log-normal distribution 
function. We choose this function because the parabola is 
the simplest function with a non-constant derivative, i.e. the 
simplest function which rises first and declines later. Since 
we are working in log-log space we had to go with its log-log 
analogue, the log-normal distribution function. The function 
peaks at M v ; r = 10 12 Mg. For the sake of simplicity we will 
refer to this kind of function as “parabola-like” throughout 
the paper. 

The ^inflow values of most data points in both panels had 
to be rescaled in the following way: the Vi n fl ow values of the 
data points within the left panel that were not at 2 = 2.46 
were scaled to this value according to the scaling relation 
presented in the right panel. Likewise the Ui n fl ow values of the 
data points within the right panel whose underlying galaxy 
bins were not at M v r = 10 12 Mg had to be rescaled to that 
value according to the scaling relation presented in the left 
panel. The colour bar axes in both panels indicate the values 
the galaxy bins used to have before rescaling. 

Both scaling relations were found in the following iter¬ 
ative process: first a square root power law was fit through 
the yet unrescaled data points in the right panel. This re¬ 
lation was then been used to rescale the data points in the 
left panel. Subsequently a log-normal distribution function 
was fit through the (scaled) data points in the left panel. Af- 
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anyway since they were at the high mass end of our collec¬ 
tion of galaxies were resolution issues play only a minor role. 
Unfortunately they were no ART or Arp galaxies available 
in this mass range. On the other hand we decided to neglect 
the MN data point at M v r = 10 11 Mq and z = 4.01 having 
a very high inflow velocity as an outlier since it has a very 
low mass, so it might be more prone to resolution effects. 

We are also aware of the fact that our choice of func¬ 
tional forms for equation 0 as well as figures [7] and [8] the 
square root power law on the one hand and the log-normal 
distribution on the other hand might seem to be very spe¬ 
cific and arbitrary but those are the log-log analogues of the 
simplest functional forms having constant and non-constant 
derivatives, the straight line and the parabola. It is surely 
possible to find other functional forms which would lead to 
a different equation 0, but due to the constraining power 
of the data points those alternative functional forms should 
ultimately arrive at very similar predictions (i.e. the almost 
same appearance of figure [8J but they will have more free 
parameters. 


6 CONCLUSIONS 


Figure 8. Average inflow velocity of the gas in units of the virial 
velocity as a function of halo mass and redshift. Shown as back¬ 
ground colour is the fitting function (equation [TJi with contour 
lines at 0.5, 0.6, 0.7, 0.8 and 0.9 u v j r . The actual values mea¬ 
sured from the simulations are given with colour coding within 
the open symbols as well as by the attached labels. ART is marked 
by squares, MN by circles and Arp by triangles. 


terwards the resulting log-normal distribution function was 
used to rescale the data points in the right panel. These steps 
have been repeated until convergence was reached (i.e. un¬ 
til the resulting scaling relations did not change anymore). 
Convoluting the square root power law with the log-normal 
distribution function leads to the following equation which 
serves as a model describing the behaviour of the cold stream 
infall velocity as a function of redshift and host halo mass: 

Uinflow(A/ v ir, z) _ A yjZ T I 

Uvir (M vil ,z) a (Mvir/M 0 ) 

x q ,|_IW.I->f| (7) 

Best fit parameter are A = (1.25±0.13) x 10 16 , = 45.4±4.8 

and a = 4.20 ± 0.53. A visualisation of this 3D functional fit 
is shown in figure [H] where the values that are actually mea¬ 
sured from the simulations are given by the colour within 
the open symbols as well as by the attached labels. The ART 
results are represented by squares, the Arp results by tri¬ 
angles and the MN results by circles. Shown as background 
colour is the fitting function (equation 0 with contour lines 
at 0.5, 0.6, 0.7, 0.8 and 0.9 iVir- 

We are aware of the fact that the choice of a “parabola¬ 
like” function instead of another power law as a scaling rela¬ 
tion for the data points in the left hand side panel of figure 
01iinges only on the four very high mass (M ^ 5 x 10 12 Mq) 
data points of the MN simulation which is the one having 
the lowest resolution. We decided to trust these data points 


In this paper we looked at the radial inflow velocities of 
accretion along streams from the cosmic web into massive 
galaxies at high redshifts using three sets of AMR hydro- 
cosmological simulations. We calculated free-fall profiles as 
well as sound speeds of the hot ambient medium via two 
different methods and we found the following: 

• The velocity profiles for the cold streams in the simu¬ 
lations are very different from free-fall. 

• Not even the infall velocities of the dark matter parti¬ 
cles follow free-fall, since there is no meaningfu l dark matter 
accretion within ~ 2 r v j r o f the central halo (|Cuesta et al.l 
120081 : IWetzel &: Nagaill2014T ) . Therefore the dark matter par¬ 
ticles are not inflowing they are rather moving back and 
forth. 

• The sound speed of the ambient medium is a better 
proxy for the behaviour of the gas inflow velocity, at least 
at high redshifts. We present a neat scaling relation for the 
inflow velocity based on this as a first order approximation 
(^inflow — 0.9 Uvir)- 

• The velocity profiles for the cold streams are by and 
large constant with radius or only very slowly increasing 
with decreasing radius. 

• There is a dependence of the shape of the gas infall ve¬ 
locity profiles with mass and redshift: for low mass haloes 
(M v i r ^ 5 x 10 12 Mq) at higher redshift (z it 1.7) one usu¬ 
ally sees an increase of velocity with decreasing radius down 
to ~ 0.5 r v i r from where the velocity is decreasing again. 
This originates in the l oss of angular moment um during the 
“strong-torque phase” dDanovich et al.ll2015l l. 

• For very high mass haloes (M v i r > 2.5 x 10 12 Mq) the 
inflow velocity is slightly increasing with increasing radius 
over the whole radius range. In all other cases the inflow 
velocity is constant or very mildly increasing with decreasing 
radius also over the whole radius range. 

• This constant infall velocity has in units of the virial 
velocity as a function of radius a “parabola-like” dependency 
on the host halo mass that peaks at M v i r = 10 12 Mq and it 
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also follows a square root power law relations with respect 
to the redshift (^inflow oc V 'z + 1 u v ir). 

• Since the magnitude of the mass accretion rate is also 
constant with radius one can follow that cold and hot gas 
do not mix as the cold streams fall in. 

A potential limitation of our simulations could arise 
from the artificial pressure floor used in the simulations 
to ensure that the Jeans mass is always fully resolved. 
It might affect the temperature and density of the very 
dense and cold parts of the streams, with implications on 
the computed inflow velocities. Also the interaction be¬ 
tween outflows and inflows is yet to be studied in simu- 
lations with strong feedback. I ndeed preliminary results of 
I House fe Dekel | in preparation! - ) show that there is no strong 
effect of strong feedback on to inflows. But AMR codes are 
still the best available tool for recovering the stream prop¬ 
erties. With their 17 pc resolution, and with proper cooling 
below 10 4 K, these simulations provide the most reliable de¬ 
scription of the cold streams that are available so far. 

We conclude that the velocity profile of the gas flowing 
into a galaxy’s halo in the form of cold streams is, contrary 
to what might be expected, roughly constant with radius in¬ 
stead of free-falling. The potential energy of the gas which is 
lost on its way is not converted into kinetic energy but must 
b e dissipated by ot her mechanisms, such as Lya radiation 
(iGoerdt et al.ll20icl l. 
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